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ABSTRACT 


The modal attenuation constants J n a cylindrical waveguide coated with a 
lossy dielectric material are studied as functions of frequency^ dielectric 
constant, and thickness of the dielectric layer. A dielectric material best 
suited for a large attenuation is suggested. Using Kirchhoff’s approximation, 
we also studied the field attenuation in a coated waveguide, which is illumi- 
nated by a normally incident plane wave. For a circular guide which has a 
diameter of 2 wavelengths and is coated with a thin lossy dielectric layer 

( " 9.1 - j2.3, thickness ■ 3% of the radius), a 3 dB attenuation is achieved 
within 16 diameters. . 
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1 . INTRODUCTION 


Reducing the radar cross section (RCS) is one of the major problems in 
designing a modern military aircraft. When an airplane is heading toward the 
radar site, a major contribution to RCS comes from the jet engine intake. The 
RCS from the jet intake is mainly due to the rim diffraction and interior irra- 
diation. Th^ rim diffraction has been stalled by several authors [IJ, [2]. 

The main goal of our research is to reduce, as much as possible, the 
interior irradiation from the jet Intake. One way to achieve this goal is to 
coat the interior wall of the jet Intake with a lossy dielectric material. Once 
the wave is transmitted from the outside illumination, the wave will attenuate 
as it propagates through the Interior of the jet Intake before it scatters back 
to the outside of the jet intake. 

For our theoretical model, we approximate the jet Intake by a cylindrical 
waveguide. We will investigate the properties of the wave attenuation in a 
cylindrical waveguide coated with a lossy dielectric material and suggest how 
the power attenuation of the transmitted wave to the waveguide from the outside 
Illumination can be maximized. 

This report begins with the derivation of the normal modes in the lossy 
waveguide. It is followed by the general discussion of the behavior of 
the attenuation constant as functions of the frequency, the dielectric constant 
and the layer thickness of the dielectric material. A few specific materials 
are chosen to show how the wave attenuates within the waveguide from the nor- 
mally Incident plane wave. In the conclusion and discussion section, other 
possible devices for a large power attenuation of the wave are suggested. 
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2. FORMULATION 


Consider a cylindrical waveguide coated with a lossy material as shown in 
Figure 1. We assume that region I is free space and the permeability of region 
II is the same as that in free space. Ui the past, a number of authors treated 
the problem of the partially filled waveguide 13], [4J. In this report, we 
rederlve the formulation to make this report self-sufficient and uniform in 
notation for other derivations presented later. 

2.1. Propagation Constant 

2.1.1. Approximate solution 

Though the perturbation theory does not give a very accurate result for the 
waveguide perturbed by a very lossy material, this analytic result provides 
guidance in the exact numerical calculation. 

The difference between the propagation constants of the perturbed and 
unperturbed waveguides is given by convention,) [ 5] 

0) L dS 

k - k -2-:: (2.1) 

® / (f X it - f X it ) «z dS 

•'S o o' 

Here ^(it) and are the fields of the perturbed and unperturbed waveguides, 

respectively, m is the angular frequency, and Ap and Ae are the differences of 
the permeabilities and permittivities between the perturbed and unperturbed 
cases, respectively. The integration is over the cross-sectional area of the 
waveguide. In this report, we assume that Ap * 0 and As * ^ ” ^^^o 

is the free-space permittivity and is the dielectric constant of the lossy 
dielectric. 
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In the cylindrically symmetric geometry, mode coupling is largest between 
the TE and TM modes with the same mode indices, e^g, , TE,j and Though the 

normal mode in this case is no longer TE or TM, it is closer to one of the two 
modes when the thickness of the dielectric layer is small. We call this mode 
"quasi** TE or TM mode and will use the same notation as in the unperturbed 
waveguide for convenience. Some authors prefer to use and EH instead of TE 
and TM, respectively[4J . 

Using the static approximation, we obtain (Appendix 1) 
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Here u Is the free-space permeability; and are the n*-*' zeros of the 
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Bessel function of order m, J and its derivative J , respectively: k^ is the 

' m m’ ' po 

radial wave vector for the unperturbed case; and 
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+ mJ^(x) 


(2.4) 


2.1.2. Exact solution 

The characteristic equation for the propagation constant of the normal mode 
in a lossy hollow cylinder can be derived by imposing the boundary condition on 
the perfectly conducting surface and matching the fields between regions I and 
II (Figure 1). The characteristic equation to be solved numerically for the 
propagation constant k^ is given by (Appendix 2) 
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Figure 1. A cylindrical waveguide coated with a dielectric material 
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liifjre and £2 are the permittivities; ^2 permeabilities; and 

k^ 2 > radial wave vectors of regions I and II, respectively; and a and b are 

the radii of the air region and the conducting cylinder, respectively. is 

the Bessel function and N is the Neumann function of order m. 

m 

2.2. Fields of the Normal Modes in a Lossy Waveguide 

Once we find the eigenvalues for the propagation constants, the eigenvec- 
tors for the fields in the lossy waveguide naturally follow. The electric and 
magnetic fields are given by (Appendix 2) 
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Here the superscipts I and II Indicate regions I and II (Figure 1), and 
subscripts p» (j) and z indicate the radial, angular and propagation-directional 
components of the fields, respectively. A, B, C and D are the constants, which 
are determined by the boundary conditions and the normalization condition. 

Those constants are related by 
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There is no mode coupling between the TE and TM modes for m * 0. We note that 
there are two degenerate taodes for each angular mode index m except when m « 0. 
In the expressions of the fields, we have arbitrarily chosen one of those 
two modes. 

2.3. Wave Propagation into a Cylindrical Waveguide from the Incident Plane Wave 

Consider a plane wave incident on the opening of a cylindrical waveguide 
(Figure 2). An exact solution for the cylindrical waveguide without the 
dielectric coating has been derived by Weinstein [6j using the Wiener-Hopf 
Method. GTD has been applied by several authors [2J, [7]. When the waveguide 
is coated with lossy material, the problem is much more complicated. Since we 
emphasize the wave attenuation within the waveguide, we use Kirchhoff's approxi- 
mation. This method provides an approximate solution and is much simpler than 
the GTD or Wiener-Hopf Method, 

We assume in the following derivation that the perturbation within the 
waveguide is weak enough so that the modal fields in the perturbed waveguide can 
be approximated by the modal fields in the unperturbed waveguide. 

For the incident fields, we write (Figure 2) 

®o “ ^ ®o^ expl-j(k^x + (2.22) 
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where 0^ is the Incident angle, k^^ and k^ are the components of the wave vector, 
is a constant, and is the free-space Impedance /y^/ e^» Note that we 
cboose the coordinate for the incideni: field in order to simplify the for- 
mulation. When the tangential E field at the opening of the waveguide is 
matched, the transmitted transverse fieldis at z » 0 are given by (Appendix 3) 
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Here k is the free-space wave number 0)/u e , and if and \f Indicate the nor- 

malized TE and TM fields, respectively. The superscripts V and H indicate the 

vertical and horizontal polarizations, respectively. The symbols with an over 

bar indicate TM modes and these without an over bar, TE modes. N and TT are 

mn mn 

the normalization constants (Appendix 3). 

When we match the tangential magnetic field, the expressions for the 
transmitted transverse fields at z 0 are similar to the above expressions 
except for the cos 0^ v.erms in Eqs. (2.28) and (2.29) and the factor of 
k /k (Appendix 3). In this report, we use the electric-field matching. 

In order to characterize the pOwer attenuation of the transmitted wave, we 
approximate the modal fields in the lossy waveguide by the modal fields in the 
perfect waveguide with the exception of the z-dependence of the propagation 
constant, which characterizes the wave attenuation. Then the power propagating 
within the lossy waveguide from a normally incident plane wave is approximately 
given by (Appendix 3) 


P(z) 
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where 
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Here is the power incident on the area of the opening of the waveguide and 
and are the real and imaginary parts of the propagation constant, 
respectively. 

We note in Eq. (2.30) that the higher mode usually carries less power 

I 

because the higher mode has a smaller but a larger 
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3. NUMERICAL RESULTS 


Numerical results are given for the dominant and modes for the 

purpose of comparison. The cutoff frequencies of lower-order modes in terms of 
the cutoff frequency of the dominant mode are shown in Figure 3. The mode pat- 
terns of 30 lowest modes are shown in Figures 4 and 5. 

3.1. Mode Attenuation and Dispersion 
3*1«1* Frequency dependence 

Figures 6 and 7 show the real and imaginary parts of the propagation 
constants of the TEj^^ and TM^^ modes as a function of frequency in a waveguide 
coated with a thin dielectric material. In Figure 6, the exact numerical solu- 
tions for the attenuation constants are compared with the results obtained from 
the perturbation theory (Section 2.1.1). We can see that the perturbation 
theory is valid only at the low-frequency region even though the thickness of 
the dielectric layer is small. At the high-frequency region, the TE^^ mode 
shows much hii'iiar attenuation than the TM^^ mode. This is due to the fact that 
the TEj^j modal field moves closer to the surface of the waveguide than the 
TMj^jl^ mode as frequency increases. These features are shown in Figures 8 and 9, 
where the magnitudes of the angular and radial electric fields are plotted as a 
function of radial distance. Note that the fields of both TE^^j^ and TMj^j modes 
at the low frequency (3 GHz) are similar to those for the unperturbed case, but 
the modal fields of the TE^^ mode (Figure 8) are closer to the surface than 
those of the TM^^ mode (Figure 9) at the high frequency (7 GHz). On the other 
hand, the real part of the propagation constant is not much different from that 
for the unperturbed case (Figure 6). 

3.1.2. Dielectric constant dependence - thin layer 

Figures 10 and 11 show the attenuation constants as a function of the 
complex dielectric constant of the lossy dielectric material. We observe two 
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Figure 3. 



The ratios of cutoff frequencies of lower-order modes to the cutoff 
frequency of the dominant mode (TEii) in a circular waveguide. 
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Figure 6. Attenuation constant as a function of the frequency. 
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Figure 8. Radial and angular coraponentc: of the electric field of the mode 

relative to the fields at p = 0: 5:2 “ “j2» a=* 9#7 cm, 

b = 10 cm, a) f * 3 GHz, b) f = 7 GHz. 
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Figure 9. Radial and angular components of the electric field of the TMi i 

model relative to the fields at p=0; £ 2 '- 1.5 - j2, a = 9.7 cm, 

b = 10 cm, a) f = 3 GHz, b) f » 7 GHz. 
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intsresting l;eatures in these graphs. First, there is a clear "resonance" 

effect of the imaginary part of e^, on the attenuation constant when the real 
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part of of is small (< 1.5). Second, a smaller gives a larger at ten- 

p 

nation constant except for that with a "dielectric resonance." Since the 

of the practical materials are usually larger than 1.5, we need to choose the 

R I 

lossy dielectric material with a large loss tangent (small but large e^) for 
a large power attenuation in the waveguide. 

Note that the general dependence of the attenuation constant on the 
dielectric constant for both and are very similar, and these proper- 

ties may not be limited to those two particular modes. To understand these 
« 

results, consider the following. The attenuation constant is proportional to 
the power dissipation within the lossy dielectric layer. The power dissipation 
per unit length is related by 


T, 'X) /■ 

“ 0 i 


dS 


e-layer 

where the integration is over the dielectric region, 
constant is given by 

/^laverl 


Im k 


2 f i X it* «z dS 


Then the attenuation 


(3.1) 


where the integration in the denominator is over the cross-sectional area. 
Hence for a large attenuation constant, we need a large and large field con- 
centration within the dielectric layer. These properties can be illustrated 
from the field distributions within the waveguide (Figures 12 and 13). The 
tangential electric fields are usually small in the dielectric region with a 
thin dielectric layer becau'.y the tangential electric fields vanish on a 
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Figure 12. Radial and angular components of the electric field of the TEj^]^ mode 
relative to the fields at p » 0: f = 3 GHz, a) £j- 1.5 - jO.4, 

b) Ej. = 1.5 - j2.2, c) Ej. = 1.5 - j8. 




perfectly conducting surface. These features are shovm In Figures 12 and 13 as 
well as the field expressions in Section 2.1.2. The ratios of the tangential 
fields to the radial fields are approximately given by 


“ x/b (3.2) 

"p 



for TM(EH) 


(3.3) 


“ ^o2*^ T TE(HE), m t 0 (3.4) 

p ^ 

Since ^p 2 X and x/b are small for a thin dielectric layer, these ratios are 
usually small. Thus E is responsible for most of the power loss. Since 
Ep at the dielectric region is inversely proportional to |e^j (see eq. (A1.7)), 
increasing decreases the electric-field strength within the dielectric layer 
(Figures 12 and 13). Since the power loss is also proportional to | 1 1 there 

may exist an optimum value of for a maximum power loss and a maximum atten- 
uation constant (Figures 12b and 13b). 

We expect the attenuation constant to become vanishingly small as becomes 
very large, because the dielectric layer becomes a perfect conductor in this 
limit. These features are shown in Figures 14 and 15, where the attenuation 
constants in the asymptotic limit of are plotted as a function of e^. Note 
that the attenuation constants decay rather slowly as becomes very large. 
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Figure 14. Asymptotic behavior of the propagation constant of with the 

variation of the imaginary part of ej.( e^. = 1.5 - jx). 








3.1.3. Dielectric-thickness dependence 


When the dielectric layer is thick, the behavior of the propagation 
constant is different from that in the case of the thin dielectric layer. These 
features are shown in Figures 16, 17, 18, and 19. We observe a few interesting 
points in these figures. 

First, there may exist a '*spatial"-resot,ance effect as the layer 
thickness Increases. That is, there may exist an optimum layer thickness for a 
large attenuation constant. As shown in Figure 20, the optimum thickness 
results when the thickness of the dielectric layer is about 1/4 of the wave- 
length within the dielectric layer. This is contrary to the common belief that 
the thicker the lossy layer Is, the larger the attenuation constant becomes. It 
is interesting to note that when the resonance is weak (e.g., TE with a small 
e^), the field is similar to the field of a surface mode confined to the 
dielectric region near the surface. This feature is shown in Figure 21a, where 
the ratioD of the magnitudes of the electric fields to those at the center of 

the waveguide are plotted as a function of the radial distance. 

I R 

Second, when we keep the ratio of e to e (i.e., the loss tangent) 

r r 

constant, the basic dependence of the attenuation constant on the dielectric- 
layer thickness remains similar. Generally, a thicker layer attenuates the TE 
mode more significantly than the TM mode. 

Third, as increases, the attenuation constant of the TM^^ mode 
Increases, but the attenuation constant of the TEj^^ mode decreases. This can be 
understood from the field distributions as shown in Figures 21 and 22. Note 
that when increases, the electric field of the TE^j mode at the dielectric 
region decreases significantly while the electric field of the TMj^j^ mode remains 
relatively unchanged. Since the attenuation constant is also proportional 
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OUTER RADIUS, b(cm) 


Figure 16, Attenuation 


constant with the variation of b(TEn,e5 



ZO*-» — i3>C3m— 4— (I> 



OUTER RADIUS, b(cra) 


Figure 17, Attenuation constant with the variation of b(TEn,£§ 
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the variation of b (TMn,e§ » 3). 





Figure 19. Attenuation constant with the variation of b (TMjj[,e^ 1.5). 
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1/4 WAVELENGTH-DISTANCE LINE 


(a) 





(c) 


(d) 



OUTER RADIUS, b(cm) 



Figure 20. The phase distances of the normal modes in wavelength in the radial 
direction within the dielectric layer (phase = Real(kQ2)(b"’a)/ r) 
with the variation of b(a = 9.7 cm, f = 3' GHz): a) TEj^^, = 3 

(Figure 16), b) TE^ , e| = 1.5 (Figure 17) „ c) TM|i , = 3 (Figure 

18), d) TMii, = 1.5 Figure (19). 
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RADIAL DISTANCE (cm) 


RADIAL DISTANCE (cm) 


Figure 21. Radial and angular components of the electric field relative to 

the fields at p * 0 (thicker layer case): a = 9.7 cm, b = 14 cm, 
a) TE^, Cj. == 1 .5 - j , b) TE^i, % * 1.5 - j2^ 
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RADIAL DISTANCE (cm) 


RADIAL DISTANCE (cm) 


Figure 22. Radial and angular components of the electric field relative to the 
fields at p 0 (thicker layer case); a = 9.7 cm, b = 14 cm, 
a.) TMii, Gr = - j, b) TM^, € 4 . = 1.5 - j2. 
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Ej. (Equation (3.1)), the attenuation constant of the mode Increases 

while the attenuation constant of the mode decreases as the e^ of the lossy 

dielectric increases. 

3.2. Wave Attenuation iu the Lossy Waveguide from a Normally Incident 
Plane Wave 

As shown in Section 3.1.2, the lossy dielectric with a large loss tangent 
is a good choice for the coating material for a large wave attenuation within 
the waveguide. Plastics are in this category [8j; three materials are chosen 
for further analysis (Table 1). Figures 23, 24 and 25 show the power attenu- 
ations of the transmitted waves from the normally incident plane wave with a 
unit power on the aperture. Only two modes (TE^^ and TEj^ 2 ) ®re propagating in 
this particular geometry (Figures 2 and 3) and 84% of the incident power on the 
aperture is transmitted in this approximation. There are two interesting 
features to be observed. Most of the power is carried by the dominant mode 
(TEii)(whlch has been discussed in Section 2.3) and the attenuation constant of 
the dominant mode is usually larger than that of the higher mode. As those two 
modes propagate through the waveguide, eventually the higher mode will carry 
most of the power, but when this happens, the total power of the wave has 
already decayed to a small fraction of the initial transmitted power at z = 0. 
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TABLE 1 


THE DIELECTRIC CONSTANTS OF THE LOSSY DIELECTRIC MATERIALS AND THE PROPAGATION 
CONSTANTS OF TE^ AND TE^2 WHEN THESE t-IATERIALS ARE USED IN THE WAVEGUIDE 

(a ■ 9.7 cm, b » 10 cm, f » 3 GHz (b/X "1)) 


Figure 

23 

24 

20 

Material 

Polystyrene 7% 

Catalin 
700 base 

Pyralin 

*Real 

9.1 

4.74 

3.74 

*Imag 

-2.275 

-0.7252 

-0.6171 

Real k xb 
TEj 

6.1321 

6.1059 

6.0960 

Imag k xb 

-1.0388x10"^ 

-6.1856x10“^ 

-7.2436x10"^ 

Real k xb 

3.3726 

3.3554 

3.3513 

Imag k xb 

TEi^ 

-9.2734x10"^ 

-2.9211x10“^ 

-2.6319x10"^ 


^Reference (8) 




































Figure 24. Attenuation along the waveguide for a normally Incident plane wave 
(catalln 700 base ( e =® 4,74 — jO.73), a = 9.7 cm, b = 10 cm, 
f = 3 GHz). ^ 







4. CONCLUSION AND DISCUSSION 


We have calculated the attenuation constants of the normal modes in the 
waveguide coated with a lossy dielectric material. When the lossy-dielectric 
layer is thin, choosing the dielectric material with a large loss tangent 
results in a large attenuation constant of the normal mode (TE or TM). We have 
chosen a few practical materials for the lossy dielectric to demonstrate the 
power attenuation of the wave in the lossy waveguide from the normally incident 
plane wave. Unfortunately, the power attenuation may not be sufficient in a 
practical application. For a circular waveguide coated with the thin layer 
((b - a)/b ■ 3%) of the best dielectric material available (polystyrene 70% and 
carbon 30%, - 9.1 - j2.3), a 3 dB attenuation can be achieved around a 

distance of 16 diameters. 

We have also calculated the attenuation constant as a function of the 
thickness of the lossy dielectric layer. There may exist an optimum thickness 
of the lossy dielectric layer for large attenuation, and the behavior of the 
attenuation constants of the TE and TM modes as a function of the dielectric 
constant of the lossy material are different from that in the thin-layer case. 
That is, choosing the lossy material with a smaller results in a larger 
attenuation constant for the TE mode but a smaller attenuation constant for the 
TM mode, or vice versa. Usually, TE modes attenuate more than TM modes in the 
thicker layer case. 

The main reason we can not obtain a large attenuation constant in the wave' 
guide coated with a thin lossy dielectric layer is that the electric field |^j 
is small because its tangential component vanishes at the perfectly conducting 
surface and its normal component is inversely proportional to the dielectric 
constant, which is usually large for the available materials. A thick lossy 
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dielectric layer may be used for a large attenuation constant of the normal 

mode. For example, when the waveguide Is coated with the thick dielectric layer 

((b - a)/b ■ 10%) of Catalin (700 base, e « 4.7 - jO.7), a 3 dB attenuation can 

r 

be obtained within a distance of one diameter. However, too thick a layer may 
not be desirable in the design of the structure. 

It may be possible to attain a large attenuation constant of the normal 
mode even with a thin dielectric layer but with a different pattern of coating. 
Consider a waveguide which is coated with double layers (Figure 26a). If we use 
the dielectric with a large dielectric constant for the outer layer, then the 
modal field will shift to the surface, and the electric field in the inner 
dielectric layer with a large will be large, where the power of the wave is 
dissipated. We have seen in Section 3.1.2 that the tangential field does not 
contribute much to the power dissipation in the case of a single thin dielectric 
layer. Increasing the tangential component of the electric field with a multi- 
layer coating makes this component play a major role in the power dissipation of 
the guided wave. The above effect with the multilayer can be achieved with 
commercially available resistive sheets (Figure 26b). 

Another way by which the electric field in the dielectric region can be 
increased may be using the corrugated layers (Figure 27a). The basic idea in 
this device is to increase the radial component of the electric field in the 
dielectric region by making space v;ithin the dielectric layer. 

Since the electric field at the center region is usually large, it has been 
suggested that putting the lossy material (e.g., resistive cards) at the center 
region may increase the attenuation constant (Figure 27b). 

The feasibility of those devices mentioned above also depends on the actual 
design problem of the jet intake. 
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(b) 


Figure 27. Other possible devices for a large attenuation constant: 
a) corrugated layer, b) resistive card. 
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APPENDIX 1 






OF POOR 


APPROXIMATE SOLUTION OF THE PROPAGATION CONSTANT IN THE 
CYLINDRICAL WAVEGUIDE COATED WITH A LOSSY DIELECTRIC MATERIAL 


In this appendix, we derive the approximate solution based on the pertur- 
bation theory. Refsr to the main text for the notations. 

The fields in the perfect cylindrical waveguide are given by 

A « zij; (magnetic vector potential) 


‘E 


' .. 1 3 ^ 

p j me 3 p3z 


“p •pT? 


f 13^ 

j mep 3<j)3z 


% 3p 


'z j me 


.2 ^ “ 


H « 0 
z 


for TM 


(Al.l) 


and 


F » z4< (electric vector potential) 


1 3ij» 


$ 3p 


E = 0 
z 


H 


H. = 


H. 


1 

2 

jmu 

3p3z 

1 

2 

jmpp 3(})3z 

1 

(7 

il- + 

j OJp 

3z^ 


for TE 


(A1.2) 
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The wave functions satisfying the boundary condition are 





f 


COS m(|) 

‘^m 

-r p 

» 4 


sin mi|) 


f > 


- 



COS mif) 

— p 

k 4 


sin m<t> 




(A1.3) 


(A1.4) 


where and are n^^ zeros of dj^(x) and J^(x), respectively. The disper- 
sion relation is 


2 

k 


z 


'h 


2 


k 


P 





where 


k 

P 


’mn 


b 


or 


t 



(A1.5) 


(A1.6) 


We approximate the perturbed fields in the denominator of Eq. (2.1) by the 
unperturbed fields. In the numerator, we use the quasistatic approximation [5] 
such that 


g = ^ I (A1.7) 

e o 
r 

where is the complex dielectric constant of the coating material. This may 
be a good approximation assuming the tangential electric field is small near the 
surface and the electric field is nearly normal to the surface in the dielectric 
region. Assuming Ap = 0, Eq. (2.1) becomes in this approximation 


k - k 
z zo 


% /d.r. 


o o 


2 / (^ X ) »z dS 
S ° ° 


(A1.8) 
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OF t'tfc.* V 

where the Integral in the numerator is over the cross section in the dielectric 
region (d.r,), and the one in the denominator is over the cross-sectional area 
of the cylinder. 

Substituting the fields foi. the TE mode in (A1.2), we obtain 


k “ k 
z zo 


L_^!2 for HE (TE ) 

Te^R. _ . mn mn' 


"r"zo 5 

/ ““ g (X) xdx 


(Al.9) 


where 


X 


(Al.lO) 


Using the recurrence formula [9], 


Substituting Eq. (Al.U) in Eq. (Al.9), we obtain 


ij? u e^( - 1) 


k - k 
z zo 


o o r 
^ ®r*^zo 


1 - 


f 


• 

i. 

b 


,2, - J 

2 


"m^ WJ 


mn' mn"^ 


(Al.U) 


( 2 . 2 ) 


where f(x) is given in Eq. (2.4). For the TM mode, we follow a similar proce- 
dure with the fields given in Eqs. (Al.l) and (A1.3). Then we obtain 


k - k 
z zo 


ZO r 


1 - 


5 a 
. mn 




I- for EH^„ (TM^„) 
' mn mn 


(2.3) 
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ORiGfwa. rv ■ ■ -s 
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EXACT NUMERICAL SOLUTION OF THB PROPAGATION CONSTANT 
IN THE CYLINDRICAL WAVEGUIDE COATED WITH A LOSSY DIELECTRIC MATERIAL [10] 

The field expressions In terms of the wave functions are given in Appendix 
1* The wave functions in Region I (Figure 1) satisfying the boundary condition 
at the origin are given by 

® (A2‘a) 

T ""jk 2# 

® ^ (A2.2) 

and the wave functions in Region II satisfying the boundary condition at the 
conducting surface ?re given by 

lV^(kp2P) - C[J^(kp2P) N„<kp2'’) " N„<‘‘p2P> ■^m^‘"p2'’^l ® 

**m^‘‘p2‘’^ " V'p2''> V'‘p2'=>* ® 

where 

2 2 9 2 

“pi + “z ’ “ ” “l 

“p2 “z * “^^>'2 “ “2 

tl/'s with and without the bar indicate TM and TE modes, respectively. For other 
notations, please refer to the main text (Section 2.1). Note that the angular 
terras of the wave functions are chosen such that these two modes are coupled. 
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Substituting Eqs. (A2.1) through (A2.4) In Eqs. (Al.l) and (A1.2), we obtain the 
general field expressions given In Eqs. (2.13) through (2.18). Matching the 
tangential fields between Regions I and II, we have four equations, 


E : 
z 


. 2 „ , 2 „ 

k ,F, k „Fo 

A_£Li. c-^ 

1 h 


“z = 


2 2 
^ iFi k „F< 

B -ELi - D-^ 


H,: Ak ,F, +-rrA-F, - Ck .Fq + — =- F 




pl^^l ^ "1 '-*"p2"3 ‘4 


E.; 


Ak_m 


. F, + Bk ,F, r, f UK^ryt/, 

(j) ojEj^a 1 pi 1 ^ p2 4 


» Ck^m 


F, + Dk.oF/ 


(A2.6> 


where the notations are given In Section 2.1.2. For a nontrivial solution for 
A, B, C and D, the determinant for Eq. (A2.6) must vanish. This condition gives 
the characteristic equation given In Eq. (2.5). 
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APPENDIX 3 




«■ 4 4 ■ 


NORMAL MODES PROPAGATING INTO A CYLINDRICAL 
WAVEGUIDE FROM THE INCIDENT PLANE WAVE 


The geometry for this problem is shovm in Figure 2, and the notations are 
A;.iicated in the main text (Section 2*3). 

The transverse electric and magnetic fields in the waveguide at z - 0^ are 
given by 




z - O"^) " I 
m.n 


mn mn mn mn k mn mn k mn mn 

0 o 


(A3.1) 


*' 4- 1 n 

- 0 ) I 


^C'' (z xf ) +^C» (z xS ) 
k ran mn k mn 

o o 


o m^n 

+ Cz xf ) (z“ x^« ) 

fMM ' MW ' mw • fnw ' 


mn 


mn 


mn mn 


mn 


(A3. 2) 


where C's are constants. The symbols with the bar indicate the TM mode, while 
those without the bar indicate the TE mode, and if and if are given in Eqs. (2.26) 
and (2.27). We can also write if and if such that 


mn mn 


z V. ’ 
t ^ran 


(A3. 3) 


if'^’^a -N 


mn 


, U) 
mn t ^mn 


(A3. 4) 


where 


^mn 


' 1 
^n 

"cos m(|)' 

^mn 

“ ‘^m 

— P 

sin m>}) 


■- 7 V 1 
%n 

^mn 



^ 'I 

^n 


sin m(j) 

" ‘^m 

— P 


cos m<}) 


. . 



(A3. 5) 


(A3. 6) 
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Define the orthonormalization condition such that 


» z 6 ,6 , 

mn m’n’ o mm’ nn' 


(A3. 7a) 


> » z 5 ,6 , 

mn m’n’ o mm' nr,' 


(A3. 7b) 




-mn -n..n-- - (A3.7c) 

mn m ' n 

Superscripts V, H indicate that these conditions apply to both cases of vertical 
and horizontal polarizations. However these two normal fields are orthogonal. 
Substituting Eq. (A3. 3) in the right-hand side of Eq. (A3. 7a), we obtain 

+V,H . 2 . " „ ,V,H. " „ V,H .. 

- <7 /’“. 7 ,/:“,> 

mn t ^mn t ^m ' n ' 


„2 fr „ ,V,H „ V,H 
N / f V*. ij^ , f dS 

mn t%n t%'n’ 


= // V? dS + i /■» i|i2- (A3.8) 

mn ' 'mn t '^ni'n' ” Imn 3p 

S c 

Due to the boundary condition at the surface, the second terra vanishes. Since 


y.H 


(A3. 8) becoraes 


^’n’ 


,V,H 
’^m'n' ' 


mn mn 



( ' 




^ , ' 


2 

^n 

2 

b 


^n 


N 6 ,6 . (1 + 6 n) TT 
mn mm' nn' ' mO' 

“ir! 
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Q. 

O 

J 

m 

b 

^ J 

P 


TT ,2, 


^2 


Then, the normalization constant becoraes 
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Equations (A3. 7b) and (A3. 7c) can be similarly proved, 
constant of the TM mode is shovm to be 

First » we assume the tangential electric field is continuous at z *■ 0: 

^£<z - o'*") (A3. 11) 


i^"(z - O') 


From Eq. (2.22) , 


1“^ (A3.9) 

The normalization 


^^”(z “ 0 ) = E^ X cos 6^ exp[“jk^xj (A3. 12) 

and the right-hand side of Eq. (A3. 11) is given in Eq. (A3.1). Multiplying both 
sides of Eq. (A3. 12) by the integration operator, 


b 2 tt „ 

/ pdp / d* if ; • 
o o 


(A3. 13) 


and using the orthonormalization conditions of Eq. (3.7), we obtain 


y 

mn 

H 

mn 


b 2 tt 

Eq cos / pdp / d^ 


0 


0 




(A3. 14) 


where 


^(x) = 


" -jk X 


X e 


Araitay and Galindo evaluated the integral of Eq. (A3. 14) (Eq. (6) in reference 
[llj). The coefficients in Eq. (A3. 14) then become 
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(A3. 15) 
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^mn 
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cos e N (“j)““^ J (5I ) J (k b) 

Z 0 mn k„ m'^n' m' x ' 

'0 

-1 

^mn 

- 

0 X 


We can following a similar procedure for the other two coefficients. They are 
shown to be 
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mn 

Vo . _ 

1 

1 

“ 0 * 


“ 7TT 

Z k_ 0 mn 

f* 'i 

- K 


G 

mn_ 

0 zmn . 

k^ - 


2 X m ^n' m X 
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(A3. 16) 


If we match the tangential magnetic field, 


H^^(z - O") = itj.( 


z ■ 0^) 


(A3. 17) 


where 

*in, ^o " -jk Jc 

H (z»0)=-^ye ^ 

and the right-hand side of Eq. (A3. 17) is given in Eq. (A3. 2). If we follow a 

similar procedure as before, the result is different from the one obtained by 

matching the tangential electric fields. The only difference in this case is 

that the cos 9 terms in Eqs. (A3. 15) and (A3. 16) and the factor of k /k . 

o zmn o 

The coefficients in this approach is given by 


mn 

.H 

^mn 


o o 

^o^zmn 
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mn k 


2tt m j ( / ) j (k b) 
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-1 


(A3. 18) 
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(A3. 19) 


In this report » we use the electrical-field matching. Witt and Price [12J indi- 
cated that using the magnetic-field matching gives a better convergency for the 
high (evanescent) modes. However, as shown in Eq. (A. 18), this method gives a 
very large coefficient for the mode near the cutoff frequency. In this case, 
the total power transmitted may be more than 100%, which can not be justified in 
a physical point of view. 

When 0Q ■ 0 (normal incidence), matching the electric fields at z » 0, 
the transmitted tangential electric field is given by 
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2En I f « — 0 
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(A3. 20) 


Note that we have terms with m » 1 only because the incident field is linearly 

A 

polarized in the x-direction. Also, from Eq. (A3, 2), the transverse magnetic 
field is given by 
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(A3. 21) 


Assuming that the perturbation is not too large, we use the fields in the 
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unperturbed waveguide for the fields In the perturbed waveguide except for the 
exponentially decaying factor In the z-dlrectlon. Then the power attenuation of 
the transmitted wave from the normally Incident field Is given by 


o"b\ “ 
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X / pdp 
0 


21‘^1 


f t 


In 


f » 




’In 


2 I 




2 , 


(A3.22) 


where 


‘‘zln * ®ln - J“ln 

P ■■ 

0 0 o 

Here Is the power Incident on the area of the aperture. The Integration In 
Eq. (A3. 22) Is evaluated In Eq. (A1.9). 
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